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Given the abstract wave equation — A a = 0, where A a is the Laplace 
operator with a point interaction of strength a, we define and study W a , 
the associated wave generator in the phase space of finite energy states. We 
prove the existence of the phase flow generated by W a , and describe its most 
relevant properties with particular emphasis on the associated symplectic 
structure and scattering theory. 
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I. INTRODUCTION. 



To introduce the problem we begin with a well known example. Given the 
free scalar, zero mass, wave equation 

4>-A<p = o, (1) 

the usual attitude in the literature is to search the solutions in the real 
Sobolev-Hilbert space H 2 (R 3 ); in order to fix the notations we recall that 
H S (R 3 ), s G R, is defined as the set of tempered distributions with a Fourier 
transform which is square integrable w.r.t. the measure with density (1 + 
\k\ 2 ) s . This is a standard mathematical choice but not the more natural one. 
In fact, equation (1) can be written in the first order form 

■ip = Wt/j , (2) 

where the linear operator 

W : H 2 (R 3 ) © H\R 3 ) -> H\R 3 ) © L 2 (R 3 ) 

is defined as 




Here A is the usual Laplace operator viewed as a self- adjoint operator on 
L 2 (R 3 ). It is well known that equation (2) generates a strongly continuous 
one parameter group of evolution 

U* : H\R 3 ) © L 2 (R 3 ) -> H\R 3 ) © L 2 (R 3 ) . 

This group group is energy preserving, i.e. there exists an energy form 

mi) = \{U\\i + \\y^<l>\\l) , 

coinciding with the Hamiltonian of the system, preserved by the flow. More- 
over U* constitutes a group of canonical transformations w.r.t. the symplectic 
form 

^ (0> 4>), 0)) : = (0> 0) - (<p, 4>) 

((•, •) denoting the usual scalar product on L 2 (R 3 )) and W is nothing but the 
Hamiltonian vector field corresponding, via cu, to £. 
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As the form of Hamiltonian S suggests, a more natural domain for the 
study of the system described by (3) is the space of the finite energy states, 
which is larger than the original one, because the first component <fi of such a 
state is not necessarily square integrable, as instead is implicit in the standard 
Sobolev environment recalled above. This more suitable description goes as 
follows. 

Let us define (general and more complete definitions will be given in the 
following section) iJ 1 (R 3 ) as the completion of the space C£°(R 3 ) in the norm 
|| \J — A 0|| 2 . Now it is possible to define the new operator W on ff^R 3 ) © 
L 2 (R 3 ), the Hilbert space of finite energy states, by 

W : # 2 (R 3 ) ®H\R 3 ) -> H\R 3 ) ©L 2 (R 3 ), W{<j>, j>) := (0, A0) . (4) 

where 

# 2 (R 3 ) := {0 G H l (R 3 ) : A0 G L 2 (R 3 )} . 

It is an easy matter to verify that W is a skew-adjoint operator (see e.g. [1, 
thm. 2.1.2], [2, §XI.10]) so that due to Stone theorem it defines a strongly 
continuous one parameter group of evolution 

U l : H\R 3 ) © L 2 (R 3 ) -> H\R 3 ) © L 2 (R 3 ) 

which is trivially energy preserving, just because the energy coincides with 
the norm of the Hilbert space, and the flow is given by a group of isometric 
operators. This procedure generalizes to the case in which one considers an 
abstract wave equation with a positive self- adjoint operator in the place of 
-A (see [3], [4, §8]). 

Here we consider and study in detail the case in which —A is replaced 
by — A Q , the Laplace operator with a point interaction of strength a (see 
section II for its precise definition), and construct the corresponding wave 
generator W a ; since — A a is not positive when a < one can not directly 
use the results appearing in [3] and [4]. 

The abstract wave equation corresponding to A Q , i.e. 

0-A Q = O, (5) 

was introduced for the first time in [5]. There, when is vector- valued and 
when a = — ™ (m the phenomenological mass, c the velocity of light, e the 
electric charge), it is shown that (5) describes the evolution of the electromag- 
netic field self-interacting with a point particle in dipole approximation (the 
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so called linearized Pauli-Fierz model). Another model connected with the 
wave equation (5), often studied in the fifties' and sixties' literature on exact 
models in quantum field theory, is the so called "pair theory" (see ([6] -[8] 
and references therein). The classical version of this model is the regularized 
version of the one we study here, and many at the time unanswered questions 
about its behaviour in the ultraviolet limit find their rigorous collocation in 
the present work. 

In [5], [9], [10] it is also shown that the Cauchy problem is well posed on 
the phase space D 1 (R 3 ) © L 2 (R 3 ), D 1 (R 3 ) ~ H 1 ^ 3 ) © R, (refer to section II 
for the definition of D 1 (R 3 )) and that the corresponding strongly continuous 
one parameter group of evolution 



where F a denotes the bilinear form corresponding to the self-adjoint operator 
— A Q . Therefore, analogously to the case of the free wave equation, the 
problem of defining (5) on the larger space of finite energy states naturally 
arises. The theory of delta point interactions was originally developped in 
the context of nonrelativistic quantum mechanics (see [11] and references 
therein); this made natural to use L 2 (R 3 ) as the underlying Hilbert space 
and so, in order to define the dynamics on the space of finite energy states, 
one has to modify the original definition of — A a , to allow the elements of 
its domain being not square integrable. This is done in section III where 
we also show (thm. 3.1) that the operators W a here constructed generate 
an evolution group a fact that, in the case a < 0, is not immediately 
evident. So, as an aside result, a conserved energy form exists; this form 
however is not positive when a < 0, and therefore it is not suitable to define 
the norm of the appropriate phase space. 

In section IV we treat the Hamiltonian formulation of the wave equa- 
tions with delta interactions. Here we solve the problem by giving a complex 
structure J a commuting with the operator W a . This leads to an equivalent 
Schrodinger-like first-order formulation which, also in view of a future quan- 
tization of the dynamical system under study, plays a key role. The complex 
structure before mentioned is obtained considering separately the case a < 



Ua '■ -D 1 (R 3 ) 



© 



L 2 (R 3 ) -> D\R 3 ) ©L 2 (R 3 ) 



preserves the energy 
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from the other case: in the strictly negative case we obtain an invariant 
splitting of the phase space and the complex structure in such a way that 
the Hamiltonian vector field appears separately as a Schrodinger equation 
both on the stable and unstable part of the phase space; in particular, on 
the unstable subspace, which is finite dimensional, the Hamiltonian is that 
of an harmonic repulsor, and the Schrodinger equation is the corresponding 
ordinary differential equation, as expected. 

The last topic treated (see section V) is the scattering theory for the pair 
of operators (W a , W). Being the Hilbert phase spaces for W a and W different, 
resort has to be made to the two Hilbert space scattering theory introduced by 
Kato in [4]. Using then Birman invariance principle and the trace condition 
of Birman-Kuroda theorem, we are able to prove the existence of the Moller 
wave operators and their completeness (thms. 5.1 and 5.2). As a consequence 
of the machinery needed for the definition of wave operators, one obtains 
a relation (see (12)-(14)) between the evolution group (acting on the real 
Hilbert space of states with finite energy) generated by W a and the unitary 
group (acting on the complex Hilbert space Lq(R 3 ), the complexification of 
L 2 (R 3 )) generated by \J — A a (in the case a < one consider only the positive 
part of the operator). This can be seen as a variation of the procedure applied 
in section III in the case one uses the standard complex structure on L^(R 3 ): 
indeed the two structures are related in a simple way (see (10) and (11)). The 
relations (12)-(14) could also be used to define the group [7*, the generator 
of which is easily seen to be W a , so providing an alternative proof of the 
existence of the dynamics. 

II. PRELIMINARIES. 

We start by giving definitions and main properties of the Sobolev type spaces 
needed in the sequel, and to which we made reference in the introduction. 
We define the family of pre-Hilbert spaces H S (R 3 ), s G R, as the set of 
tempered distributions with a Fourier transform (denoted by A or by T) 
which is square integrable w.r.t. to the measure with density \k\ 2s . The 
scalar product is defined as 



Note that, when s > 0, H S (R 3 ) C H S (R 3 ) and H' S (R 3 ) C H' S (R 3 ), the 
embeddings being continuous. Since \k\~ 2s is locally integrable for any s < 
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3/2, 

2' 



Vs<-, L 2 (R 3 ,\x\ 2s dx) c S'(R 3 ), H S (R 3 ) = ^-\L 2 (R 3 , \k\ 2s dk)) 



and thus H S (R 3 ) is complete for any s < 3/2 and coincides with the usual 
Riesz potential spaces (see e.g. [12, §7.1.2]). 

We can then define the isomorphism (r — s < 3/2) 

(-A) 3 / 2 : iT(R 3 ) -> iT-*(R 3 ) , ^((-A) s / 2 0)(A;) := |A;| s 0(A;) . 

Our notation is justified by observing that, in the case < r = s < 3/2, 
(-A) s / 2 coincides with the closure of (-A) 5 / 2 : H S (R 3 ) -> L 2 (R 3 ). 

Since, contrarily to what happens for the usual Sobolev chain if s (R 3 ), 
^T*(R 3 ) is not included in H s (R 3 ) when r > s, we also define the sequence of 
spaces 

n / n \ 

#«(R 3 ) : = p| H k (R 3 ) =^- l [ p| L 2 (R 3 , |x| 2fc dx) . 

fc=i \fe=i / 

Obvioulsy H n (R 3 ) is a Hilbert space with norm 

/ n \ V 2 

||0bn:=^Ell(-A) fe/2 0lliJ • 

We come now to point interactions; for their general theory of we refer to 
[11]; here we confine ourselves to the essential definitions and results. The 
operator — A a describing a standard point interaction at the origin with 
strength a is defined as follows. Let us introduce the dense linear subspace 
ofL 2 (R 3 ) 



G L 2 (R 3 ) : = A + Q^ A , 0aG# 2 (R 3 ), (a + ^Q^MO)), 



where < A ^ — sign(o;) (Ana) 2 and 

-\f\\x\ 

The Laplacian with a point interaction with strength a is the operator 
-A a : D 2 a (R 3 ) - L 2 (R 3 ) , -A a( f> := -A0 A - \Q^G X . 
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Its resolvent is given by 

(-A a + A)- 1 = (-A + Ar 1 +^ + ^) 



where G A ® G A (0) := (G A , 0) G A . 

The bilinear form corresponding to — A Q has domain D 1 (R 3 ) x D 1 (R 3 ), 
D l {R 3 ) := {0 G L 2 (R 3 ) : = A + Q G A , A G /^(R 3 ), G R} , 
and is defined by 

F a (0, <p) := ((-A + A) 1/2 A , (-A + A) x /Va) - A(0, ip)+(a + ^ 

(see [13]). Both the expressions for F a and — A Q contain the arbitrary param- 
eter A, but contrarily to the appearance, they do not depend on it. Indeed 
(following [14, §2]) the operator and form domain can be defined in the fol- 
lowing alternative, and more useful, way. Note that, since for any A > 

G\ G L 2 (R 3 ) , G-G A Gtf 2 (R 3 ), (G - G x ) (0) = ^ , 

where 

G(x) [ 



+ Q<t>{G\ — G) G H 2 (R 3 ) 



4:7l\x\ ' 

defining 

0reg •" 

we have equivalent ly 

^(R 3 ) = 

[<P G L 2 (R 3 ) : = <p reg + Q^G, <j) reg G H 2 {R 3 ), Q G R, aQ^ = re9 (O)} . 

Correspondingly, the form domain is 

D\R 3 ) = {0 G L 2 (R 3 ) : = <j> reg + Q^G, re9 G tf^R 3 ), Q G R} 

so that, with this definition, the singular part of the field is exactly Coulom- 
bian. However such a singular field G is not in the configuration space D 1 (R 3 ). 
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The removal of this incongruence will lead, in the following section, to the 
introduction of the operator W a . 

With the domains so given we can redefine the operator and the form as 

-A a (f) = -A(f) reg 

and 

F a (<j),<p) = ((-A) 1/2 (f) reg , (-A) 1/2 ip reg ) +aQ (j) Q (f . 

Now it is well known (see [11, Chap. 1.1]) that — A a is a selfadjoint operator 
in L 2 (R 3 ). An important property is that — A a is positive only for a > 0, 
whereas for a < it is only bounded from below; more precisely if a > 
(repulsive delta interactions) the spectrum of the operator is absolutely con- 
tinuous and coinciding with [0, +oo); if a < (attractive delta interactions) 
the spectrum is given by {— A } U [0, +oo), where — A = —(Ana) 2 is an 
isolated negative eigenvalue, and the remaining part of the spectrum is ab- 
solutely continuous. In the Schrodinger case this eigenvalue corresponds to 
a bound state, while in the wave case where one has a second order equation 
in time, it leads to unstable solutions exponentially running away in the past 
or in the future (see [5], [9], [10] and reference therein for the meaning of 
these well known runaway solutions in classical electrodynamics). 

We now come to the wave generator associated to the standard delta 
operator. Its domain and action are given by 

W a : Z^(R 3 ) © D\R 3 ) -> D\R 3 ) © L 2 (R 3 ) , 

*($)-U !)($)■ 

By considering the Hilbert space structure given by D 1 (R 3 ) © L 2 (R 3 ) ~ 
-f/^R 3 ) © R © L 2 (R 3 ) this operator is the generator of a strongly continu- 
ous group of operators 

[/* : D 1 (R 3 ) © L 2 (R 3 ) -> D 1 (R 3 ) © L 2 (R 3 ) . 

In the case a > this is an immediate conseguence of the skew-adjointness 
of W a with respect to the positive energy scalar product on the phase space 
given by 

(((0, 0), {tp, tp))) a := (0, 0) + F a (0, ip) . (7) 

More precisely one has the following result (the proof being a straighforward 
calculation): 
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Theorem 2.1. For any a G R, with respect to the scalar product ((•,-)) / g, 
/3 > 0, one has 

D(W* a ) = D(W P ) 

and 

W*{4, 0) = - (j> reg + ^ Q^G, , (3 > , 

W*M, 0) = ~ (0re 9 + Q^G, A0 re9 ) = W^(0, 0) , /3 = = a . 



In the case a < the operator W Q is readily proven to be a generator by 
considering the operator 

W a , A (0,0) := (0, (A a - A)0 re9 ) , 

where A > A . This, being now — A Q + A positive, is skew-adjoint with respect 
to the scalar product 

(((0,0),(^)» QjA := (0,0> + F Q (0,v?) + A(0,^> , (8) 

and so it generates a group of isometries (w.r.t. the Hilbert structure given 
by (8)). The original operator W a , being a perturbation of the previous 
one by a bounded operator, also generates a strongly continuous group of 
operators on the phase space (which however are no more isometries). 

We now describe an alternative way to prove that W a , a < 0, is a gen- 
erator. Such a different method will play a key role in the next sections. 
As we already said before in the case a < the self-adjoint operator — A a 
has a negative eigenvalue — A (with corresponding normalized eigenvector 
47T\/— 2a G\ ) which gives rise to the runaway solutions of the wave equation 
associated to W a . Proceeding as in [10, §4] (note that there we worked with 
the different decomposition = (f>\ + Q<pG\ ) we consider the linear operator 

w? : [£«(R 3 )W © [-D 1 (R 3 )]nr -> PWW © [£ 2 (R 3 )W , 




where, given any vector subspace V C L 2 (R 3 ), we have defined the corre- 
sponding "non runaway" subspace [V] nr by 

[V] nr :={0GV : (0,G Ao ) = O}, 



9 



and 

A" r := (A a )|[ £) 2( R 3)] TH . = P nr ■ (A a )|[ D 2 (R3)]nr , 

P nr being the orthogonal projector onto [L 2 (R 3 )] nr . By simple calculations 
one has (see [10, §4]) 

[D 1 (R 3 )] nr = j<rj G D 1 (R 3 ) : Q^ = -47r V ^(0 reff ,G Ao )} 
[^(R 3 )]n. = {0 G ^(R 3 ) : re9 (O) = A (0 re9! C Ao )} 

and 

K r <P 

= A<f) reg - 871-^3^ (A<p reg , G Xq )G Xo 

= A(f) reg + 8n^J\~ ((—A + A o )0 re9 , G Xo )G x „ - 8tv^J\~ \ (0 re9 , G Xo )G Xo 
= A<p reg . 

The non-negative bilinear form associated to — A" r is then 

F?{4>,<P) = ({-AY l2 <t>re 9 ,{-^) 1/2 Vre 9 ) ~ ^ (<j> r e a , G Xo ) ( Vreg , G X() ) 

and W™ r is skew-adjoint w.r.t. the scalar product 

(((0,0),(y^)>>- := <0,0> + iC(0,y.). 

The strongly continuous one parameter group of evolution generated by W™ r 
preserves the non-negative energy 

-=\ (<0,0> + ^W)) 

with coincides with the Hamiltonian of the sistem w.r.t. the symplectic form 
u (see [10, thm. 4.2] for an alternative Hamiltonian picture). 
Since A a G X() = \ G X() , and 

D 1 (R 3 ) ~ [D 1 (R 3 )] nr © R , 
L 2 (R 3 )^[L 2 (R 3 )] nr ©R, 
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we can write 

W a = W2 r xAo, 

where 

A :R 2 ^R 2 , A (x,x) := (x,\ x) . 
Therefore W a , a < 0, is a generator and 

U* a = e tw ° = e tw « x e tA ° . 

Here and below, given two linear operators Ai : D{A{) — > Hi and A 2 : 
D(A 2 ) -> i? 2 , x i 2 : x D(A 2 ) ^Hi®H 2 denotes the the linear 

operator defined by 

A 1 x A 2 (0!,0 2 ) := (Ai0i, A 2 2 ) ■ 

In conclusion, for any a G R, £7* is a group of canonical transformation w.r.t. 
the symplectic form a;, and W a is the Hamiltonian vector field corresponding 
to the energy 

£ a (j>, 4) = l - ( + F a (0, 0) ) = £(0 re „ 4) + 1 Ql . 

Let us remark that the flow [/* can be explicitly calculated (see [5, thm. 
3.1]). 

III. THE OPERATOR W a . 

Now we would like to mimic the construction of the energy space for the 
usual wave generator and the extension of the operator itself, to the case of 
delta point interactions. To this end, let us define the linear operator 

W a : Dl (R 3 ) © D 1 (R 3 ) - D 1 (R 3 ) © L 2 (R 3 ) , VT a (0, 0) := (0, A0 re9 ) , 

where 

D 2 a (R 3 ) := {0 = re , + Q^G, <p reg G H 2 (R 3 ), Q<f, e R, aQ (t> = re9 (O)} , 

D\R 3 ) := {0 = re9 + Q^G, re9 G ^(R 3 ), g G R} . 

Analogously to the free case /^(R 3 ) © L 2 (R 3 ) describes now the space of 
finite energy states. Moreover the Coulombian singularity G is now in the 
configuration space Z) 1 (R 3 ) . 
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Introducing the Hilbert space structure given by Z) 1 (R 3 ) © L 2 (R 3 ) ~ 
H 1 (R 3 ) © R © L 2 (R 3 ) we want now to show that also in this case W a gener- 
ates a strongly continuous one parameter group of evolution. When a > 0, 
considering, similarly to the case of W a , the scalar product 

(((0, 0), (ip, tp))) a := (0, y?> + ((-A) 1/2 re9 , (-A) 1/2 ip reg ) + a Q^Q V , 

one can prove that W a is skew-adjoint and so it is a generator. Note that 
when a — 0, contrarily to situation discussed in the previous section, ((•, -)) a 
is no more a scalar product, being annihilated by the zero energy eigenvector 
(G, 0) (this fact has to be compared with the presence of a zero energy 
resonance for — A ). In order to show that also in the case a < W a 
is a generator one can not use the same strategy as before consisting in a 
translation, since the scalar product (8) is now ill-defined, I) 1 (R 3 ) being not 
a subset of L 2 (R 3 ). So the perturbation argument fails and we are forced to 
proceed in an alternative way. The decomposition of W a , a < 0, introduced 
at the end of the previous section is our starting point: we simply extend it 
to the case of W a . Therefore we define, when a < 0, 

[D 2 a (R 3 )] nr = {<!> e D 2 a (R 3 ) : re9 (O) = A o (0 re9 ,G A() )} 
[D\R 3 )] nr = {0 G D 1 (R 3 ) : = -AiT^X, (0 reg ,G Ao )} , 

and 

W2 r : [Dl{R 3 )U © [D\R 3 )] nr - [D\R 3 )] nr © [L 2 (R 3 )] nr , 
W™ := (0,Ar0), 

where 

A-:[L> 2 (R 3 )] nr ^[L 2 (R 3 )] nr , 
A^ r 0:= A0 re9 . 

With such definitions W™ r results skew-adjoint with respect to the scalar 
product on [D 1 (R 3 )] nr © [L 2 (R 3 )]„ r given by 

{{{M)M<P)))a ■= (0,^) + ((-A) 1/2 re9 ,(-A) 1 /V e9 ) 

-47T \l /2 ((j>reg,G\ )((p re g,Gx ) ■ 

Moreover, since 

D 2 a (R 3 ) ~ [ J D 2 (R 3 )] nr ©R, D\R 3 ) ~ [£) 1 (R 3 )]„ r © R, 
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similarly to the case of W a we have 



W a = W^xA . 



For the case a = 0a similar decomposition is possible by using the projection 
onto the subspace orthogonal to the eigenvector (G, 0). Indeed, defining 



W {0) : Hq(R 3 ) © -D^R 3 ) - ^(R 3 ) ©L 2 (R 3 ) , W {Q) {4>,4>) := (0 re9 , A0) , 



where #o(R 3 ) := \0 e # 2 (R 3 ) : 0(0) = Ok the operator W(q) is skew-adjoint 



(((0, 0), (ip, 0))) (o) := (0, 0) + ((~A) 1/2 0, ("A) 1 ^) 
and, since Z) 1 (R 3 ) ~ /J 1 (R 3 ) © R, the following decomposition holds: 

W = W {0 ) x . 



We can now state our result regarding the existence of dynamics: 

Theorem 3.1. W a is a closed operator coinciding with the closure ofW a . 
It generates a strongly continuous group of evolution 



where {0n}i° C D 1 (R 3 ) is any sequence such that n — > in D 1 (R 3 ). 

Proof. W a is a generator since it is skew-adjoint when a > and Wq = 
W(o) x 0, W a = W™ r x A , a < 0, where both W( ) and W^ r are skew- 
adjoint. Therefore W a is closed. By its definition W a is equal to W a on 
D 2 a (R?)®D l {R?) and so it coincides with the closure of W a if Dl(R 3 )®D 1 (R 3 ) 
is a core. This is proven as follows: 

analogously to the case of W a , any e Z) 2 (R 3 ) admits the representation 



with respect to the scalar product 



[7* : D\R 3 ) © L 2 (R 3 ) -> ^(R 3 ) © L 2 (R 3 ) . 



which can be defined as 



C7*(0,0) = hmC/*(0 n ,0), 



= 0A + Q<f>G\ 



where 



0A = <Preg " Q*(G A - G) G ^(R 3 ) 
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and 




where 




Now if 0^ converges in H 2 (R 3 ) to 0a, we have that Q n converges to Q<f,, thanks 
to the continuous embedding of H 2 (R 3 ) in C fe °(R 3 ) (see e.g. [12, §5.6.2]). 

Being W a equal to W a on _D 2 (R 3 ) © D 1 (R 3 ) the same is true for the cor- 
responding groups of evolution. Since L> 2 (R 3 ) (BD 1 (R 3 ) is dense in D 1 (R 3 ) © 
L 2 (R 3 ) which is continuously embedded in Z) 1 (R 3 )©L 2 (R 3 ), one has the equal- 
ity £/*(0,0) = C/*(0,0) for any G D 1 (R 3 ). The proof is then concluded by 
the denseness of D 1 (R 3 ) in /^(R 3 ). □ 

Let us remark that, since W a is the closure of W a , our construction coincides, 
in the case a > 0, with the abstract one given in [3] (see also [4, §8] for a 
similar construction). Moreover, since W a = W™ r x Ao when a < 0, one has 
that W™ r is the closure of W™ r - 

IV. THE SYMPLECTIC STRUCTURE. 

The standard symplectic structure recalled in the introduction, 



it is not well defined on the phase space finite energy states, i.e Z) X (R 3 ) © 
L 2 (R 3 ). This requires a different approach to the Hamiltonian description 
of the dynamical system described in the previous paragraph. The problem 
shows up already in the case of the free wave equation, with the phase space 
i/ 1 (R 3 ) © L 2 (R 3 ); usually in the standard literature on infinite dimensional 
Hamiltonian systems (see e.g. [15]) only the easier case of the free field with 
strictly positive mass is explicitely discussed. 

We recall that (see [16], [15]) when the Hilbert space carries a complex 
structure J7", it is possible to complexify the space in such a way that the 
imaginary part of the complex scalar product turns out to be a symplectic 
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form, while the real part is the old (real) scalar product, coinciding with 
the energy. Any skew-adjoint operator A commuting with J remains skew- 
adjoint within the complex Hilbert space, so that %A := J ■ A is self-adjoint. 
Therefore, since e tA = e~ lt ^ lA \ A generates a strongly continuous group of 
unitary (hence symplectic) transformations. More precisely, collecting the 
known results on the subject (see e.g. [15, §2.6, §2.7], [16, chap. II]), we 
state the following: 

Theorem 4.1. Let A be an injective skew-adjoint operator on the real Hilbert 
space H with inner product (■,■). Then the closure of the densely defined 
linear operator 

A ■ {-A 2 )' 1 ' 2 : Range (A) -> H 

defines a complex structure J commuting with A. Defining, for any ip G H , 
the multiplication by the complex number i as 

iifj ■= Jifj ; 

H becomes a complex Hilbert space with Hermitean inner product 

[if>i, ^2] := (^1, ^2) + i (1P1, J1P2) ■ 

The strongly continuous one parameter group U l := e tA is a group of sym- 
plectic transformations relatively to the symplectic form 

0(^i, ^2) := Imf^i,^] 

and the linear vector field 

A : D(A) -> H 

is Hamiltonian with associated densely defined Hamiltonian function 

where Q denotes the quadratic form associated to the self-adjoint operator 
J -A. 

A wide class of examples is obtained by the following construction, which 
is a simple consequence of the above theorem. Let us consider an injective 
nonnegative self-adjoint operator 

B : D (B) — > K 
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on the Hilbert space K and let us consider the closure of 

( M 

\-B 2 ) 

on the Hilbert space H = D(B) ®K, where D(B) is the completion of D(B) 
with respect to the norm \\u\\b '■= \\Bu\\k- In the case in which this closure 
is injective, the complex structure J given by the previous theorem is 

J B : D(B) D{B) © K , J B (u, v) = (B^v, Bu) , 

where B and B^ 1 are the closures respectively of B and its inverse B^ 1 : 
Range (B) — > D(B). This allows to endow H with the structure of a com- 
plex Hilbert space, which we continue to call H; precisely, defined a generic 
element as w := (u,v) G D(B) © K, the Hermitean scalar product in H is 

[w 1 ,w 2 ]b ■= ((wi, w 2 ))b + i ({wi, Jbw 2 ))b , 

where 

((w 1 ,w 2 )) B ■= (Bu!,Bu 2 ) + (vi,v 2 ) ■ 
On the product H x H we have the symplectic form 

n B : D{B) ®Kx D(B) ®K^R tl B (w 1: w 2 ) = ((w u J B w 2 )) . 

With respect to the complex variable w the wave equation 

u = -B 2 u (9) 

assumes the Schrodinger-like form 

—iw = Bw . 

Moreover such an equation is Hamiltonian w.r.t. the symplectic form Qb 
and the densely defined Hamiltonian function 

H B : D(B^ 2 ) x D(B l l 2 ) - R, H B {w) = \ (}\B l ' 2 v\\ 2 + \\B^u\\ 2 ) , 

where the operator B s is defined as the closure of B s . 

The strongly continuous symplectic group of operators obtained by solv- 
ing the equation (9) preserves the energy £ B (u,v) := \ [w,w]b- 
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An immediate example is given by the choice B = \/—A : if 1 (R 3 ) — > L 2 (R 3 ), 
corresponding to the standard wave equation and leading to the complex 
structure 

J : H\R 3 ) © L 2 (R 3 ) -> ^(R 3 ) © L 2 (R 3 ) , 

:= ((-A) _1/2 0, -(-A) 1/2 0) . 

Other concrete examples are obtained when the operator B 2 is a point inter- 
action, more precisely B = ^f^K^ : D 1 (R 3 ) -> L 2 (R 3 ) with a > 0. In this 
case the corresponding complex structure is given by 

J a : ^(R 3 ) © L 2 (R 3 ) -> L> 1 (R 3 ) © L 2 (R 3 ) , 

^(0,0) := ((-^"^^-(-Aa) 1 ^) • 

The same procedure is not directly applicable to the cases a < 0, due to the 
lack of skewadjointness and injectivity for the operator W a . A natural way 
out is to project the operator on the subspace of absolute continuity and to 
apply the abstract scheme to this projection. This works well for the case 
a < 0, whereas the case a = deserves a different treatment. Here are the 
details of the two constructions. 

In the case a < we have seen in section III that W™' is skew-adjoint, w.r.t. 
the scalar product ((•, •)) 1 a T , and one-to-one. Therefore we can apply to it 
thm. 4.1 (or better the successive example with B = ^J— A™ r ) obtaining the 
complex structure J a commuting with W a , a < 0, defined as 

Ja '■= J a x 3 j 

where 

J? : [D\R*)\ nr © [£ 2 (R 3 )W - [-D 1 (R 3 )]nr © [^ 2 (R 3 )]n, , 

Ja{4>A) = ((-KT 1/2 t,-(-K r ) 1/2 <f>) 

and 

j : R 2 -> R 2 , j(x, x) := (x, -x) . 
Here, analogously to the case a > 0, the linear operators 

(-Ar) 1/2 :[^(R 3 )W-[L 2 (R 3 )] nr , 
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and 

( _ Arr i/ 2;[L2(R3)U ^ [5 l (R 3 )U) 

are defined as the closures of 

(-A™') 1/2 : PWU C [D\R*)] nr -> [L 2 (R 3 )]„ r 

and 

(-A-)- 1 / 2 : Range ((-A" r ) 1/2 ) C [L 2 (R 3 )] nr - [5 1 (R 3 )] nr 
respectively. 

We have then the complex Hilbert space of the couples 

z) := ((0, 0), (x, x)) E [D\R 3 )] nr © [L 2 (R 3 )]„ r © R 2 , 
with the Hermitean scalar product 

[Ol, Zl), (lp 2 , Z 2 )]a ■= ifa]™ + [Zi, Z 2 \ , 

where 
and 

[z u z 2 ] := {z 1 ,z 2 )+i{z 1 ,jz 2 ), (zi,z 2 ) := x 1 x 2 + x x x 2 . 
The associated symplectic form is 

il a : [D 1 (R 3 )] nr © [L 2 (R 3 )]„ r © R 2 -> R , 

n«((^i^i),(^,%)) = ((i>uJS r i>2)) n : + {zujz 2 ) . 

With respect to the complex variables (if;, z) the wave equation corresponding 
to W a takes the Schrodinger-like form 

-ij> = (-AT) 1/2 *P T J-Ao o\ 

-i z = L z °- [ 1 J ' 

and such an equation is Hamiltonian w.r.t. the symplectic form Q a and the 
densely defined Hamiltonian function 

H a : D(Q n J) © R 2 -> R , W a ty, z) = \ Q n JW + \ (L z, z) , 

18 



where 

qsw) = \ (ii(-Ar) 1/4 0ii^ + ii(-Ar) 3/4 0ii^) 

is the quadratic form associated to the self-adjoint operator J™' ■ W" r . 

If (0,0) G D 1 (R 3 ) © L 2 (R 3 ) has the orthogonal decomposition (0,0) = 
■0 + zG, where G denotes the normalized eigenvector corresponding to Ao, 
then 

Therefore, being 

JJ^ = e tW S r x e tA 

a strongly continuous group of unitary and symplectic transformations, the 
energy is conserved by the flow. 

We come now to the case a = 0. In this case, in order to apply thm. 4.1, 
which requires injectivity, it is necessary to project onto the subspace orthog- 
onal to the eigenvector (G, 0). Being W(p) one-to-one and skew-adjoint w.r.t. 
the scalar product ((•, -))(o), one can then apply thm. 4.1 thus obtaining a 
one parameter group Ufa of symplectic transformations such that 

U* = £/$,, x 1 

and so Uq preserves the energy £ o (0, 0) = ^(0re 9 , 0)- 

Since (— Ao)^ 2 = (— A) 1 / 2 res (note that this equality holds true only in 
the case a — 0) one has 

res = (-A)- 1 / 2 • (-A o ) 1/2 

and so, when (0,0) 6 # 2 ( r3 ) © ^(R 3 ), 

J-W {0) {4>A) = ((-A) 1 / 2 0,(-A o ) 1 / 2 0) = ((-A o ) 1/2 0, ("Ao) 1/2 0) • 

Moreover J commutes with W(o) (see (15) in the next section) and so J 
coincides with the complex structure associated to W(o) by thm. 4.1. With 
respect to the complex variable ip = (0, 0) the wave equation corresponding 
to W(p) assumes the Schrodinger-like form 

-# = (-A ) 1 / 2 ^. 
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We summarize the results obtained in the following 

Theorem 4.2. For every a G R\{0} there exists a symplectic form 

tt a : D\R 3 ) © L 2 (R 3 ) x D 1 (R 3 ) © L 2 (R 3 ) -> R 

with respect to which the vector field 

W a : D 2 a (R 3 ) x D\R 3 ) -> ^(R 3 ) © L 2 (R 3 ) 

Hamiltoniam. Moreover for a < i/ie analogous result occurs for the 
reduced vector fields 



: [5* (R 3 )]„ r x [D (R 3 )] nr - [D\R 3 )] nr © [L 2 (R 3 )] 

ano? 

iy (0) : Hq(R 3 ) © £ X (R 3 ) -> ^(R 3 ) © £ 2 (R 3 ) . 

For every a G R £/ie evolution group [/* preserves the energy 



V. SCATTERING THEORY. 

The Hilbert spaces where the operators 1T Q and IT act on, respectively 
D 1 (R 3 ) © L 2 (R 3 ) and H 1 (R 3 ) © L 2 (R 3 ), are different (also as sets), and so 
one is forced to use a two Hilbert space formulation to treat scattering the- 
ory for the pair (W a ,W). We refer to the seminal paper by Kato [4] for 
the relevant constructions and results in scattering theory with two Hilbert 
spaces. Our approach will follow the lines of the construction given in [4, 
§8-9] (also see [17, §3.5]). 

From now on, given the real vector space L 2 (R 3 ), we will denote by Lq(R 3 ) 
the complex vector space 

L 2 (R 3 ) :={0 1+ *0 2 , 1? 2 GL 2 (R 3 )} . 

We begin introducing the isometries 

C : H l (R 3 ) © L 2 (R 3 ) — > L 2 C (R 3 ) , 
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C(0,0)^C o (0,0) := (-A) 1/2 - i<j>, 
C a : ^(R 3 ) © L 2 (R 3 ) ^L 2 C (R 3 ), a >0 , 
C a {4>,4>) := (-A a ) 1/2 - i<p , 
C a : [D 1 (R 3 )]„ r © [L 2 (R 3 )]„ r - [L 2 (R 3 )]„ r , a < , 
C a ((0,0)) := (-Ar) 1/2 0-^. 
These isometries lead to the following relations: 

J = C~ 1 -iC, (10) 

J a = C- 1 -iC a , a>0, J^ = C- l -iC a , a<0, (11) 
U* = C' 1 ■ e itV ^ ■ C , 
U* = C" 1 • e**^/^ -Cxi, (12) 



U t a = C- 1 -e u ^ : -C a , a>0, (13) 



C7* = C' 1 ■ e u V^ ■ C a x e tA ° , a<0. (14) 
Note that the two equalities 

J = C~ 1 -iC, U\ Q) = C' 1 ■ e^^o • C (15) 

imply, as we stated in the previous section, that J commutes with Wmy 

Moreover the relations (12)-(14) provide an alternative construction of the 
dynamics generated by W a . In fact one could use such relations as definitions 
of [/* and then check by differentiating with respect to the time parameter 
that this evolution group is generated by the operator W a . 
We introduce now the identification operators 

J a : D 1 (R 3 ) © L 2 (R 3 ) -> H\R 3 ) © L 2 (R 3 ) , 

f ((-A)~ 1 /2.(_A Q )V20 ) 0) ; for«>0 

J Q (0,0):= (0 re9 ,0) = ((-A)- 1 /2. ( _A o) i/2 ^ )? fora = 

{ ((-A)-V2.(_A-) 1/2 -n nr 0,0), fora<0, 
where U nr denotes the projection 

U nr :D 1 (R 3 )^[D\R 3 )] nr , 
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and 

J' a : H\R 3 ) © L 2 (R 3 ) -> D\R 3 ) © L 2 (R 3 ) , 

' ((-AJ-V^-A) 1 /^), fora>0 
J^(0,0):= (0,0), fora = 

{ ( (-A^)- 1 / 2 • P nr ■ (-A) 1 / 2 ^,^) , fora<0. 

We can then define the Moller wave operators 

n ± (W, W a , J a ) := s- lim C7-* • J a ■ U f a ■ P ac (W a ) , 

t^±oo 



n ± (W a , W; J' a ) := s- Jim^ I/-' ■ J' a -U\ 



where 



P ac (W a ) : D\R 3 ) © L 2 (R 3 ) - ^(R 3 ) © L 2 (R 3 ) , 
Pac{W a ){(t>A) := < 



(0, 0) , for ct > 

(0re 9 , 0) , for a = 

k (n nr 0, P„ r 0) , for a < . 



Concerning the existence of such wave operators we have the following 
Theorem 5.1. The Moller wave operators 

n ± (W, W a] J a ) := s- lim U-* ■J a -U t a - P ac (W a ) , 

t— >±oo 

fi ± (W a , j;) := s- t hrn^ • J* a • C/* 
exist, are complete and are mutually adjoint isometries, i.e. 

Range W^; •/<*) = Range W a \ J a ) = li?i(R3)© L 2 (R3) , 
Range Q+(W a: W; J' a ) = Range Sl-(W a , W; J' a ) = Range P ac (W a ) , 

fi ± (w, j a y ■ n±(w, w a , j a ) = p ac (w a ) , 
n ± (w a , W; j' a y ■ n ± (w a , w- j' a ) = i mm ^ m , 
n ± (w,w a ,j a y = n ± (w a ,w-,j' a ). 

Proof. With the above definitions one has 

n ± (w, w a , j a ) = c- 1 ■ n±(V=A, JiF a , i a ) ■ c a ■ p ac {w a ) , 
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n±(w a , w- j' a ) = c a x ■ n±(JjT a , v^; i' a ) ■ c 



where 



Q±(v /Z A, JH a ; I a ) := s- lim e" 

v t — >±oo 



itJ^A 



At\/ H a 



■ I ■ e l 



n±(y/H a ,V=A;r a ) :=s- lim e 

v t — »zboo 

Here I' a := P ac (— A a ), I a is its left inverse, and 



-it\/H a 



■I' a -e 



itV-A 



-A, 



for a > 



( -A" r , for a < 
By Birman invariance principle one has 

Ot(>/=A, yfH^;I a )=Q ± (-A,H a ;I a ) 

and 



n±(^/H a: V^K;0=n ± (H a: -A;i' a ). 

Therefore one has the identities 



s- lim e 



s- lim e ltA ■ I a ■ e 

tA 



itH a 
itA a 



and 



= n±(-A,-A a ) 



n ± {jH a ,V=E-r a ) = s- lim e-^« • I' a . e 

= s- lim P ac (-A Q ) • < 

f^±oo 



-it A 



= P ac (-A a )-Q ± (-A a ,-A) 
= fi ± (-A Q ,-A). 

In conclusion one obtains the equalities 

Q ± (W, W a , J a ) = c- 1 ■ n±(-A, -A a ) ■ C a ■ P ac (W a ) , 

Q ± (W a , W; J' a ) = C- 1 ■ n ± (-A a , -A; ) • C 

and the proof is concluded since the wave operators Q±(—A,—A a ), and 
Q±(—A a , —A) exist, are complete and are mutually adjoint isometries. This 
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is proven (see [11, appendix E]) by the Birman-Kuroda theorem being the 
resolvent difference 

(-A a + z)- 1 -(-A + z)- 1 

a rank one (hence trace class) operator. □ 

The previous theorem holds true also with the different (a-independent and 
much simpler and natural) couple of identification operators defined by 

J : ^(R 3 ) © L 2 (R 3 ) - H^R 3 ) © L 2 (R 3 ) , J(0, 0) := (0 re „ 0) , 

J' : H\R 3 ) © L 2 (R 3 ) -> D\R 3 ) © L 2 (R 3 ) , J'(0, 0) := (0, 0) . 

This is true by [4, thms. 10.3 and 10.5] since the condition 10.1 in [4] is 
verified with m — M — 1. In our situation such a condition simply reads 

V0e# 1 (R 3 )nRangeP ac (-A a ), 11^011^ = ||v / =A0||l 2 ■ 

In more detail one has the following 

Theorem 5.2. J is ([?*, ±)- equivalent to J a , i.e. 

s- lim (J a - J).Ul -P ac (W a ) = 0, 

t^±oo 

Therefore 

fi±(Vy, W a ; J) :=s- lim U'' ■ J ■ ■ P ac (W a ) 

t— >±oo 

exzst and are egtia/ to Jl±(Vy, W 7 ^; J a ). 

J' is a (E/*,±)- asymptotic left-inverse to J, i.e. 

s- lim (J'- J-l)-U t a -P ac (W a ) = 0, 

t— >±oo 

n±(^ ai W"; J') := s- lim U' 1 -J'-U 1 

t — »ioo 

exist and are equal to Q±(W, W a ; J)* = fi±(Wa, W 7 ; J^)- 
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